Two-Stage Rotational Disordering of a Molecular Crystal Surface: Ceo 
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We propose a two-stage mechanism for the rotational surface disordering phase transition of a 
molecular crystal, as realized in Ceo fullerite. Our study, based on Monte Carlo simulations, uncovers 
the existence of a new intermediate regime, between a low temperature ordered (2 x 2) state, and 
a high temperature (1 x 1) disordered phase. In the intermediate regime there is partial disorder, 
strongest for a subset of particularly frustrated surface molecules. These concepts and calculations 
provide a coherent understanding of experimental observations, with possible extension to other 
molecular crystal surfaces. 
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Surface-specific thermal disordering phenomena such 
as for instance surface melting [[[] , are of considerable in- 
terest. Conceptually, they demonstrate how a chunk of 
the metastable phase (e.g, the liquid phase within the 
solid domain of stability) can be imprisoned in full equi- 
librium by restricting it to an interface region of space. 
Practically, they provide at a temperature below the 
triple point a germ of the disordered phase that makes a 
major difference to all surface-related properties, as well 
as to the bulk transformation kinetics. 

In molecular crystals, thermal surface disordering is 
richer than in atomic solids. Because of the presence 
of rotational degrees of freedom, the bulk develops ro- 
tational and positional melting transitions that are sep- 
arate, and both of which might possess a surface coun- 
terpart. One is thus led to seek the existence of a well 
defined surface rotational melting for a molecular crys- 
tal. However until recently, very little phenomenology 
has gathered addressing that issue. 

That situation recently changed after the discovery of 
the Ceo molecule |p[ and its subsequent availability in 
macroscopic quantities || in its crystalline form, the ful- 
lerite. Here the near sphericity of the molecules favors 
rotational disordering at a relatively low temperature, 
when positional disorder and evaporation are still irrele- 
vant owing to a sufficiently large cohesion (J] . On the the- 
oretical side, the remarkable symmetry of the molecule 
leads to considerable simplifications with the help of sym- 
metry adapted functions || , such as those originally used 
for the CD4 molecules [||. 

Bulk measurements showed that Cgo-fullerite under- 
goes a first-order rotational melting transition at 260K, 
characterized by fixed orientations of the molecules in 
the space group Pa3 at low temperature transforming 
to a plastic phase with rotational isotropy, space group 
Fm3m, at high temperature [^ |9|. In the Pa3> phase, 
the fullerite structure is composed of four sublattices. 



Neutron powder diffraction studies JlfJ of this phase 
have revealed the orientational configuration of the Ceo 
molecules. The optimal orientation is such that for all ad- 
jacent couples of Ceo, an electron- rich double bond of one 
molecule faces an electron-poor pentagon of the other. In 
the Fm3m phase, the molecules are still not completely 
free to rotate as they are still correlated, and moreover 
they feel a crystal field of cubic symmetry. However, 
in this phase angular correlations are short-ranged, and 
molecules are all equivalent, as opposed to the Pa3 phase, 
with four inequivalent molecules and long-range angular 
correlations. 

The structure of a free (Iff) surface of fullerite can 
differ from that of a bulk (III) plane due to the lowering 
of the number of neighbors. The first experiments show- 
ing a specific surface disordering were high-resolution re- 
flection electron-energy-loss spectra (EELS) obtained by 
Goldoni et al jflj . This technique is sensitive to the sur- 
face layer only due to the short escape depth of the scat- 
tered electrons. The authors reported an abrupt jump 
around 230K in the energy of the Cgo triplet exciton. 
This discontinuity is an indication that the fullerite (111) 
surface undergoes the rotational disordering phase tran- 
sition at a temperature some 30K lower than the bulk. 
A similar transition temperature (230 ± 20K) has been 
observed in LEED studies of Ceo films [[l2| but the incer- 
titude in the temperature was too large to interpret it as 
a specific surface effect. 

At first sight, such a lowered surface transition temper- 
ature is reminiscent of ordinary, positional surface melt- 
ing of solids, where surface premelting precedes and ini- 
tiates the bulk transition [jjj. However, it soon appeared 
that in fullerite the surface disordering took place in the 
topmost layer only [jll 13 1, in contrast with standard sur- 
face melting, where the liquid film grows progressively 
invading the crystal. At the same time, a conceptual dif- 
ference was also noted. Positional order in a crystal is 
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forced exclusively by interparticle interaction; rotational 
ordering instead involves both intermolecular angular in- 
teractions and a cubic crystal field that would be present 
and effective on each molecule even if all the surround- 
ing molecules were perfect spheres. A surface molecule 
thus feels not only less rotational interactions, but also 
a crystal field which may generally frustrate bulk order 
jl4]]. On the other hand, due to the weak short range 
nature of intermolecular potentials, the molecules in the 
second layer, immediately below the surface are already 
very close to a bulk situation. The system, between 230K 
and 260K, can thus be regarded as an ordered solid with 
a disordered surface. 

We now underline the fact that experimental data sug- 
gest in fact not one, but two surface transitions, the sec- 
ond one taking place at even lower temperature. There 
is a less- marked discontinuity around 160K in the triplet 
exciton energy revealed by EELS It was explained as 
a consequence of the molecular rotation relaxation time 
which at this temperature becomes similar to the exciton 
lifetime. However an independent evidence for a transi- 
tion or crossover at that lower temperature is the pres- 
ence around 160K of a shoulder in the LEED intensity 
of the Pa3 diffraction spots versus temperature, absent 
in the bulk X-ray intensity |l5| . Helium atom scattering 
(HAS) diffraction studies of the (111) surface of fullerite 
by Glebov et al clearly showed a change from a (2X2) 
to a (1X1) structure taking place at 235K. But, in addi- 
tion, the behavior of the Debye- Waller Factor (DWF) of 
the (0,0) Bragg spot between 140K and 235K was com- 
pletely unexpected. Instead of decreasing monotonically 
with increasing temperature, the DWF started increas- 
ing around 140K and reached a maximum at 235K. This 
was explained as a consequence of preempted disordering 
of Ceo molecules at low coordination defect sites. This 
behavior, found to be reversible, cannot be explained if 
the surface layer were totally ordered below 235K. 

In the present work, we studied the surface orienta- 
tional ordering of fullerite by classical Monte-Carlo (MC) 
simulation. We considered a Cqq molecule as a rigid 
body, and used the potential developed by Lamoen |17| | 
to model the intermolecular interactions. That model 
considers 210 interacting centers on each fullerene: 60 
are localized on the carbon atoms, 60 are at the center 
of the single bonds, and three are distributed along each 
of the 30 double bonds to take their spatial extension 
into account. The repulsive part of the intermolecular 
potential is a sum of Born-Mayer terms between all the 
interacting centers (c) and the attractive part is made of 
R~ 6 van der Waals interactions between the atoms (a) 
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Our MC simulations were based on a standard 
Metropolis, acceptance-rejection algorithm. The 



molecules were placed in a semi-infinite fee lattice with 
a (111) surface. We considered twelve non-equivalent 
molecules, four in the first layer (the surface layer), four 
in the second layer, and four in the third layer. All other 
molecules in the system were assumed to be copies of this 
twelve-molecules basic cell, with periodic boundary con- 
ditions applied in two dimensions. In the direction of the 
bulk, the second layer was reproduced at the fourth layer, 
shifted horizontally in order to simulate the fee stacking. 
The center of mass of the molecules were held fixed to 
lattice sites, while the temperature dependent lattice pa- 
rameter of the structure was set to the experimental value 
Est]. The twelve molecular orientations were allowed to 
change randomly. We performed a complete calculation 
of the energy at each MC step. 

Earlier theoretical studies of the rotational bulk phase 
transition have shown the usefulness of a set of order pa- 
rameters based on the cubic harmonics ||. We used an 
efficient and original algorithm to generate these cubic 
harmonics. Starting from generating functions of the ir- 
reducible representations k {A\, A 2 , E, T\, T 2 ) of the 
cubic group Oh, it is possible to expand them into Tay- 
lor series and to extract from that expansion the number 
of irreducible representations for each manifold I of the 
intermolecular potential. Using the generators T (k) of 
each irreducible representation, one can construct auto- 
matically the cubic harmonics 
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where g belong to Oh and E^ k '^ is the integrity basis 
of order I for the irreducible representation k. Since the 
cubic harmonics must be homogeneous polynomials of 
the coordinates, this system can be solved easily, and one 
obtains analytical expressions for all the cubic harmonics, 
to be used with the atomic coordinate of the 60 atoms of 
the molecule. In particular, the following expression was 
calculated 
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where r™ is the position vector of atom a in the n th 
molecule in the unit cell. Using the Ceo atomic coordi- 
nates automatically introduces the icosahedral symmetry 
of the molecules. Another advantage of this method is to 
naturally allow for a future inclusion of internal molecular 
vibrations, not treated at this stage. Due to the icosahe- 
dral symmetry, hn' k ^ differ from zero for I = 6, 10, 12, 16, 
18 • • • . We used I = 10 throughout. For each molecule 
(index n), we calculated a scalar order parameter 
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where {hn^' ^ is an average over the configurations gen- 
erated by the MC algorithm. For all the irreducible rep- 
resentations of the Oh group, this parameter was traced 
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as a function of the temperature to observe the phase 
transition as it has the largest amplitude in the ordered 
phase. 

As mentioned in p0| T\ g and Tig are the two pri- 
mary order parameters. We adopted the irreducible rep- 
resentation Tig as the best indicator of the rotational 
phase transition, as the corresponding order parameter 
vanishes in the disordered phase i|. Figure 1 shows 
the temperature dependence of that order parameter for 
the twelve molecules of our system. The eight bulk-like 
molecules of layers 2 and 3 (represented by triangles), 
plus all their infinite copies, undergo a clear phase tran- 
sition at 250K. This value is affected by size, in partic- 
ular by artificial correlations built in a small system like 
the one we considered. By comparison, the mean-field 
theory transition temperature was however even smaller, 
namely 200K |l7| . However, only the dominant part of 
the multiple order parameters was originally used in the 
mean- field molecular equations of ref . |]l7f . An extension 
of the number of order parameters did raise the mean 
field transition temperature to 230K |)| . 

At low temperature, the molecules in the (111) surface 
(represented by circles in figure 1) keep the same orien- 
tation as in the bulk, corresponding to the (2 x 2) super- 
structure observed in HAS diffraction. In this configura- 
tion, shown in figure 2, one can immediately distinguish 
two kinds of molecular orientations. Three molecules in 
the surface (2 x 2) unit cell have only a two-fold axis nor- 
mal to the surface and the Pa3 threefold axis 30 degrees 
off. The fourth molecule, which we denote as BS, has 
instead the three-fold axis normal to the surface. Calcu- 
lations performed by rotating one molecule while keeping 
a fixed orientation of all the other molecules show that 
the "crystal field" angle-dependent potential felt by the 
rotating molecule is different, and shallower for the BS 
molecule than for the three others. This suggests that 
the rotational thermal behavior of the BS molecule must 
differ from that of the other surface molecules [|l4| . 

In fact with increasing temperature, the four surface 
molecules do behave differently. The BS practically be- 
comes a quasi- free rotator above around 150K whereas 
the other three molecules do so only above a much higher 
phase transition temperature, which is around 230K. 

In the low temperature phase, the surface symmetry 
is threefold due to the local threefold axis of the BS 
molecule. At 150K, we observe a dramatic reduction of 
threefold order with strong rotational disordering of the 
BS molecule. Although strictly speaking the threefold 
order of that molecule cannot vanish, it does however be- 
come so small to assimilate it to a free rotator, in which 
case the surface symmetry would become sixfold because 
the BS molecule environment is characterized by a six- 
fold axis. Thus 150K represents a crossover temperature, 
where there is no strict change of surface symmetry and 
probably no proper surface phase transition, but there 
is nonetheless a change from essential order to essential 



disorder for BS molecules only. 

Hence our system exhibits altogether four phases or 
regimes (see table I): 

(a) a fully ordered Pa3 state (effective p3 surface sym- 
metry) below 150K. (b) an intermediate regime between 
150K and 230K with one essentially disordered and three 
ordered surface molecules per unit cell (p6 effective sym- 
metry of the surface layer) over a Pa3 ordered bulk, (c) 
a disordered surface layer on an ordered PaZ bulk be- 
tween 230K and 250K. (d) a fully rotationally disordered 
Fm3m (p6m symmetry of the surface layer) above 250K. 

In the simulations, full thermal disordering of the 
molecules takes place at the topmost layer only in agree- 
ment with evidence from EELS Jill . In addition, the 
main surface transition temperature (230K) agrees very 
well with the observations by LEED, EELS and HAS. 
The intermediate BS-disordered regime between 150K 
and 230K can explain the weaker experimental anomalies 
mentioned above, in particular the feature at 160K in the 
triplet exciton energy, signaling in our interpretation an 
approximate change of surface symmetry from p3 to p6. 

We note that in LEED and HAS diffraction, this 
change is not picked up since the unit cell remains (2 x 2) 
up to 230K. We suggest instead that STM experiments 
be carried out in the appropriate temperature interval to 
test in more detail our predicted two-stage surface dis- 
ordering. Our model provides a plausible explanation 
for the behavior of the (0,0) Debye- Waller factor in HAS 
diffraction p6| , its increase between 140 and 235K due to 
the BS molecule rotating freely in that temperature inter- 
val. Due to the resulting partial disordering of the (2 x 2) 
surface cell, the intensity of some diffracted beams is re- 
duced and transferred to the specular spot |^TJ, which 
compensates the thermal attenuation of the elastic re- 
flection [fl6t . Our proposed selective thermal disordering 
of BS molecules may also provide an explanation of the 
anomalous features observed near 160K in X-ray powder 
diffraction intensities Such anomalies were also 

observed by David in high-resolution neutron scat- 
tering. If these anomalies were not reported later on, it 
is probably due to the increasing quality (and then the 
grain size) of the Cgo powder, decreasing the surface/bulk 
ratio. 

In summary, the study of fullerite may prove of 
paradigmatic importance to a much broader class of 
molecular crystals. Our study indicates that rotational 
thermal surface melting of a molecular crystal differs 
from ordinary, positional surface melting, in one impor- 
tant aspect, with several consequences. The rotational 
ordering forces comprise both crystal fields and inter- 
molecular interactions, while the former is missing in the 
positional case. At the surface, both interactions and 
crystal fields are different from the bulk. When, as in 
fullerite, the crystal field is quantitatively very impor- 
tant in the bulk, its heavy change at the surface may 
lead to an early thermal disordering taking place strictly 
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in the first layer, as in all other layers the crystal field 
is bulk-like. Moreover, crystal fields may differ so much 
between one surface molecule and another, to cause sur- 
face thermal disordering to take place in separate stages, 
the weaker molecules disordering much earlier than the 
stronger ones. 
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TABLE I. Ordering phase transitions in C6o-fullerite with 
a (111) surface. 



temperature bulk 



surface 



surf. 



sym. 



below 150K ordered 
150K-230K ordered 
230K-250K ordered 



ordered 

1 of 4 mol. 

3 of 4 mol. 
disordered 



disordered 
ordered 



above 250K disordered disordered 



p3 

P 6 

p6m 
p6m 



FIG. 1. Variation of the order parameter ijn correspond- 
ing to I = 10 and k = versus temperature. The curves cor- 
respond to the three behaviors observed in our Monte-Carlo 
simulations. The triangles are the arithmetic means of the or- 
der parameters associated to the bulk molecules. The circles 
are the arithmetic means of all the surface molecules but the 
molecule BS which has a behavior pointed out by the squares. 
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FIG. 2. Equilibrium configurations of the C6o fullerite in 
the (111) plane with space group Pa3. The unit cell is (2 x 2). 
Three of the four molecules have a two-fold axis normal to the 
plane. The fourth has a three-fold axis. This last molecule 
experiences a six-fold symmetry interaction from its neighbors 
located in the same layer. 
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